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We propose to create superposition states of over 100 Strontium atoms being in a ground state 
or metastable optical clock state, using the Kerr-type interaction due to Rydberg state dressing 
in an optical lattice. The two components of the superposition can differ by of order 300 eV in 
energy, allowing tests of energy decoherence models with greatly improved sensitivity. We take into 
account the effects of higher-order nonlinearities, spatial inhomogeneity of the interaction, decay 
from the Rydberg state, collective many-body decoherence, atomic motion, molecular formation 
and diminishing Rydberg level separation for increasing principal number. 


I. INTRODUCTION 

There are currently many efforts towards demonstrat¬ 
ing fundamental quantum effects such as superposition 
and entanglement in macroscopic systems p3fT6]. One 
relevant class of quantum states are so-called cat states, 
i.e. superposition states involving two components that 
are very different in some physical observable, such as 
position, phase or spin. Here we propose a method for 
creating such large superpositions in energy. This is rele¬ 
vant in the context of testing proposed quantum-gravity 
related energy decoherence 

Our method relies on the uniform Kerr-type interac¬ 
tion that can be generated between atoms by weak dress¬ 
ing with a Rydberg state [2Q]42^ . This can be used to 
generate cat states similarly to the optical proposal of 
Ref. [23]. Using an optical clock state in Strontium as 
one of the two atomic basis states makes it possible to 
create large and long-lived energy superposition states. 
The superposition can be verified by observing a char¬ 
acteristic revival. We analyze the effects of relevant im¬ 
perfections including higher-order nonlinearities, spatial 
inhomogeneity of the interaction, decay from the Ryd¬ 
berg state, atomic motion in the optical lattice, collec¬ 
tive many-body decoherence triggered by black-body in¬ 
duced transitions, molecular formation, and diminishing 
Rydberg level separation for increasing principal num¬ 
ber. Our scheme significantly improves the precision of 
energy decoherence detection. 

Previous related, but distinct, work includes Ref. [24] 
who briefly discussed the creation of energy superposition 
states in Strontium Bose-Einstein condensates based on 
collisional interactions. Ref. [25] proposed the creation 
of energy superposition states of light, and ref. [26] re¬ 
ported the realization of 14-ion GHZ state, with 24 eV 
energy separation, but without mentioning the energy 
superposition aspect. The present proposal promises 
much greater sensitivity to energy decoherence thanks 
to a much longer lifetime (compared to Ref. [25]) and 
to both increased size and longer lifetime (compared to 
Ref. [26]). Related work involving Rydberg states in¬ 
cludes Refs. m [28] . who performed detailed studies 
of the creation of moderate-size cat states using Ryd- 



FIG. 1. (Color online) Proposed scheme for creation of large 
energy superposition, (a) Level scheme in Strontium. The 
pseudo-spin states are the singlet ground state \g) and a long 
lived excited triplet state |e). An off-resonant laser field (Ur) 
dresses the excited state with the Rydberg level |r). This cre¬ 
ates a Kerr-type interaction between the atoms in the excited 
state. The resonant laser field (Ue) is applied for population 
rotation, (b-d) The evolution of the Husimi distribution of 
the collective spin state on the Bloch sphere. Application of 
the Kerr-type interaction splits the initial coherent spin state 
(CSS) (b) into a superposition of two CSS at opposite poles of 
the Bloch sphere (c). Applying a 7r/2 rotation along the x axis 
following the cat creation process results in a superposition 
of all atoms being in the ground or excited state. 


berg blockade. The number of atoms is limited to of or¬ 
der ten in these schemes due to competing requirements 
for the presence and absence of blockade between dif¬ 
ferent Rydberg transitions in the same ensemble. They 
also don’t use metastable optical clock states, result¬ 
ing in only small differences in energy between the two 
components. Ref. [29] briefly discussed the creation of 
moderate-size (15 atoms) GHZ type states in Strontium 
atom chains, without mentioning the energy superposi¬ 
tion aspect. Ref. [29] uses attractive Rydberg interac- 
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tions, but not the uniform Kerr-type interaction used in 
the present work. The number of atoms in Ref. [29] is 
limited by unwanted transitions to other nearby many- 
body states [30] . 

The paper is organized as follows. We begin with a 
description of our scheme in Sec. |n| In Sec. [ni| and [IV 


we quantify the effects of the main imperfections and de¬ 
coherence sources on the fidelity of final cat state. In 
Sec. El we find an estimate for size of cat states that can 
be realized with high fidelity. We then show that our 
scheme is experimentally realizable in Sec. lYU followed 
by a detailed discussion in Sec. |VII[ demonstrating that 
the effects of atomic motion, molecular formation, col¬ 
lective many-body decoherence, level mixing and BBR 
radiation induced decoherence can be suppressed. We 
conclude the paper in Sec. VIII with a discussion of the 
application of energy superposition states for the detec¬ 
tion of energy decoherence. 


II. SCHEME 


To avoid the scattering of population from the ground 
dressed state to the excited dressed state, the coupling 
term 0 = should be smaller than 

(see realization section VI for examples). 

Focusing on the ground dressed state, the effective light 
shift of the system is 


E_ = 




A2 ^ 


(4) 


Within the weak dressing regime (^ i) one can 
Taylor expand the light shift to 


E_ = 


A 

'2 




2 A2 


1 

8 A4 


N 

-0(^)^)], (5) 


which can be simplified to E_ ^ ~ with w = 

^ and xo = 2u;^A. Therefore adiabatic weak dressing 
of atoms to the Rydberg level imposes an effective Kerr- 
type Hamiltonian 


We now describe our proposal in more detail. In an 
ensemble of N ultra-cold Strontium atoms trapped in 
a 3D optical lattice m, one can consider a two-level 
system consisting of the singlet ground state \g) and a 
long-lived excited triplet state | e ), which are separated in 
energy by 1.8 eV. An interaction between the atoms can 
be induced by dressing the clock state with a strongly 
interacting Rydberg level [2Qll22] as shown in the level 
scheme of Fig. This induces a light shift (LS) on the 
atoms which depends on the Rydberg blockade. 


H = {Ne 


2 


(6) 


on the atoms within the blockade radius. The effects of 
higher order terms in the Taylor expansion are discussed 
in Sec. ImAl and Fig.j^ 


B. Generation of Cat State on the Equator of the 
Bloch Sphere 


A. Kerr-type Rydberg Dressed Interaction 


When the entire ensemble is inside the blockade radius, 
the dressing laser couples the state with no Rydberg exci¬ 
tation = (g)^ \(j)i) (where (j) G {e^g}) to a state where 
only one of the atoms in the |e) level gets excited to the 
Rydberg level |'02) = Yhi |0i---n---0Ar) with an enhanced 
Rabi frequency [32], where Ne is the number of 

atoms in the excited state. Over the Rydberg dressing 
process, the Hamiltonian can be diagonalized instanta¬ 
neously 

(I) 

where E± = ^(1 ± and 

/ cos(6»/2) -sin(6»/2) \ . 

sin(0/2) cos(0/2) j ' 


with 6 = tan Xhe Schrodinger equation ex¬ 

pressed in the dressed state basis \ip >= U\'il> > is 


( W-) 

dt \ 


E. -id/2\f \cp_)\ 
iei2 E+ )\ W+) ) ■ 


( 3 ) 


The two levels \gi) and je^) for each atom are equiv¬ 
alent to a spin 1/2 system with Pauli matrices cfx^ = 
(l^i)(o| + \ei){gi\)/2, (7^^ = i{\gi){ei\ - \ei){gi\)/2 and 
= {\ei){ei\ — \gi){gi\)/2 acting on the atom at site 
i. We define collective spin operators Si = ^ 

coherent spin state (CSS) is defined as a direct product 
of single spin states [33] 

|6>,0) = (8),^i[cos6>|^)i + sin6>e*^|e)i], (7) 

where all the spins are pointing in the same direction, and 
(j) and 0 are the angles on the (collective) Bloch sphere. 
The CSS can also be represented as [33] 


N 


| 7 ?) = \e,4>) = Y. E^C{N,N,)\N-N,), 

Ne=0 

(8) 

where g = tan(6>/2)e C{N,Ne) = 


N 

Ne 


and 


is the Dicke state of Ne excited atoms, where \N;Ne) 
is an alternative representation of the | J M) basis with 
N = 2 J and W = J + M. 
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Let us now discuss the time evolution of an initial CSS 
17^) under the Kerr-type interaction of Eq. The state 
evolves as 

N 

Ne=0 

(9) 

At the “cat creation” time Tr = — the linear term of 

^ Xo 

Eq. ^creates a phase rotation, which changes the state to 

^ _ ■ NeXO 

= |e ^ 2 w^ ^^T]). The quadratic term produces coeffi¬ 
cients of (1) and {—i) for even and odd respectively. 
The state can then be rewritten as a superposition of two 
CSS, namely 

li’iTc)) = 1 ^') + I - v')) (10) 

in analogy with Ref. [23] . Continuing the interaction for 
another Tc, one can observe the revival of the initial CSS. 
This revival can be used as proof for the successful cre¬ 
ation of a quantum superposition at Tc, since a statistical 
mixture of CSS at Tc would evolve into another mixture 
of separate peaks naEn. 



FIG. 2. (Color online) Effect of higher than second order 
nonlinearities (from the higher orders of Eq. on the fidelity 
of the cat state. The weak dressing parameter (le = has 
to be reduced for larger atom numbers N in order to keep a 
fixed fidelity Fni {Fni = 0.7 (green), 0.8 (red), 0.9 (blue) from 
top to bottom). The inset shows the Husimi Q function for an 
A = 100 cat state (a) with Fni =0.9 (corresponding to the 
black cross in the main figure), as well as the corresponding 
revival (b). The approximate revival of the initial CSS at the 
time t = 2 tc proves the existence of a quantum superposition 
at t = Tc. 


C. Creating the Energy Cat 


III. IMPERFECTIONS 


To create an energy superposition state we thus have 
to apply the following steps. Starting from the collective 
ground state we apply a 7r/2 pulse on the \e) — \g) 

transition that results in the maximum eigenstate of the 
Sx operator |?7 = 1) = as shown in Fig.fl|b). 

Since the atoms are confined to the ground states oi op¬ 
tical lattice traps, the position-dependent phase factors 
associated with laser excitation of the clock state are 
constant over the course of the experiment and can be 
absorbed into the definition of the ato mic ba sis states 
(detailed discussion can be found in Sec. VII A). We now 
apply the Kerr-type interaction. The large coefficient of 
the linear term in the Hamiltonian leads to a rotation 
of the created cat state on the equator of Bloch sphere. 
Applying accurate interaction timing, the state can be 
chosen to be a superposition of two CSS pointing to op¬ 
posite directions along the y axis on the Bloch sphere 
\i’i.Tc)) = = i) + e-*T| 7 ; = -i)), see Fig. j^c) 

and inset (a) of Eig.[^ For example, a timing precision of 
dTc = A results in an adequate phase uncertainty 

of 5(j) = (examples can be found in the realization 
Sec. VI). Applying another | pulse on the created cat 


state results in 




-, which is a superposition of 


all the atoms being in the ground and excited states, as 
shown in Fig. [^d). The created state is a superposition 
of two components with very different energies. To verify 
the creation of the energy cat state one needs to rotate 
the state back to the equator and detect the revival of 
the initial CSS under the Kerr-type interaction, see also 
the inset of Fig. |^b). 


In this section we quantify the effects of the most im¬ 
portant imperfections with direct impact on the achiev¬ 
able cat size. Other sources of imperfections, which can 
be made to have relatively benign effects on our scheme, 
are discussed in Sec. Ivnl 


A. Higher Order Non-linearities 

First, we only considered the linear and quadratic 
terms in in our Hamiltonian, which is accurate for 
very weak dressing. Applying stronger dressing fields 
yields a stronger interaction, but also increases the im¬ 
portance of higher order terms in Eq. (§. To quantify 
the effects of these higher orders, we calculate the fidelity 
of the cat state (|'0'(rc))) generated based on Eq. Q with 
respect to the closest ideal cat state, 

Fni = max(,,0,a,Tel(^'('rc)|;^(|6',<(') +e*“|7r-6»,()) + 7r))|^. 



has to be reduced for larger atom numbers in order to 
achieve a desired fidelity. 


B. Effects of Interaction Inhomogeneities 

We also considered a uniform blockade over the entire 
medium, leading to a homogeneous interaction. In prac¬ 
tice the interaction is not perfectly homogeneous. One 
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FIG. 3. (Color online) Effect of interaction inhomogene¬ 
ity. (a) Plateau-type interaction between each pair of atoms 
dressed to the Rydberg state. The interaction is uniform for 
separations up to of order the blockade radius, (b) Infidelity 
caused by interaction inhomogeneity as a function of cat size 
(N), for a constant blockade radius. Non-linear fidelity is set 
to Fni = 0.9, the blockade radius Rb = 3.6/im is created by 
Rydberg dressing to n = 80, and the atoms are considered to 
be in a cubic trap with space diagonal D and lattice spacing 
of 200nm. 

can apply fourth order perturbation theory to find the in¬ 
teraction of the entire weakly dressed system as [35l |36] 

q2 

H = Y.x{rii)KeHe-j^F. ( 12 ) 

i<j 

The many-body interaction is the sum of binary interac¬ 
tions 

x{nj)=Xo (13) 

'ij ^ 

where Rb = is the blockade radius in the weak 

dressing regime. This binary interaction has a plateau 
type nature, see Fig. [^a). The inhomogeneity of the in¬ 
teraction introduces a coupling to non-symmetric states, 
since the Hamiltonian no longer commutes with the to¬ 
tal spin operator ^ 0). We evaluate the fi¬ 

delity of a cat state created by the realistic non-uniform 
interaction with respect to the ideal cat state. Writ¬ 
ing the pair interactions xi'^ij) terms of small fluc¬ 
tuations Cij around a mean value Xm^ we decompose 
the Hamiltonian into a sum of two commuting terms, 

Vh=E XmKeHe - sN = Xm{^^) - « 

Kj 

and ViH = y Corresponding 

i<j 

to the homogeneous and inhomogeneous parts respec¬ 
tively. While the homogeneous part leads to an ideal 
cat state, the inhomogeneous part reduces the fidelity 
by a factor Fjh = |(^ = = l)p, where 

|?7 = 1) = is the initial CSS. Taylor expand¬ 

ing the inhomogeneous part of the evolution operator one 
obtains an estimate for the fidelity as explained in Ap¬ 
pendix 1^ Fig. [^b) shows the resulting infidelity as a 
function of cat size for constant blockade radius. 


IV. DECOHERENCE 

The main source of decoherence in our system is de¬ 
population of the Rydberg level which also determines 
the lifetime of the dressed state (rg ~ In this 

section we identify different Rydberg decay channels and 
discuss their effects on the fidelity of the cat state. Loss 
due to collisions is reduced by the use of an optical lattice 
trap with a single atom per site. Ref. m implemented 
a Strontium optical clock using a blue-detuned lattice 
(trap laser wavelength 390 nm) with a collision-limited 
lifetime of 100s, demonstrating that loss due to the trap 
laser can be made negligible. Other sources of decoher¬ 
ence including blackbody radiation induced transitions, 
collective many-body decoherence and molecular forma¬ 
tion will be discussed in Sec. ivm 


A. Rydberg Decay Channels 

The main source of decoherence in our system is de¬ 
population of the Rydberg level which also determines 
the lifetime of the dressed state (rg ~ The Ry¬ 

dberg state depopulation rate can be calculated as the 
sum of spontaneous transition probabilities to the lower 
states (given by Einstein A-coefficients) [38H4Q] 

0^2 

F = = (14) 

/ Ef<Ei 

where ujif = is the transition frequency and {i\f]f) 

is the dipole matrix element between initial and final 
states (see Appendix IS)- The summation is only over 
the states |/) with lower energies compared to the initial 
state. Using a cryogenic environment [41], black-body ra¬ 
diation induced transitions are negligible, see Sec. jVHCj 
for detailed discussion. 

Considering the dressing to bsns^Si in our proposal, 
the possible destinations of dipole transitions are limited 
to ^Po,i, 2 : due to the selection rules. Around 55% of the 
transferred population will be trapped within the long- 
lived ^P 2 states, which we refer to as qubit loss. Around 
35% of the population is transferred to ^Pi states, which 
mainly decay to the ground state \g) = within a 

short time (e.g. r^sbp^Pi = 23/i5 [42|), which we refer to 
as de-excitation. The remaining 10% of the population 
is transferred to ^Pq states. Half of this population (5% 
of the total) contributes to qubit loss, bringing the total 
loss to 60 %, while the other half (also 5% of the total) is 
transferred to the excited state, which effectively causes 
dephasing of |e) because the photon that is emitted in 
the process contains which-path information about the 
qubit state. 
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B. Effects of Rydberg; Decoherence on the Cat 
State 

The three decoherence types discussed in the previous 
sub-section have different effects on the cat state. Loss 
and de-excitation completely destroy the cat state if they 
occur, while dephasing is both unlikely and relatively be¬ 
nign. We now explain these statements in more detail. 

The majority (60%) of the dressed state’s decay goes 
to non-qubit states |e) ^ (5|e)|0)p + \/l — S‘^\l)\l)p^ where 
and \ l)p represents the emitted photon. In 
addition to loss, 35% of the dressed state’s decay is de¬ 
excitation |e) ^ 6\e)\0)p + Vl — where = 

g-0.357eTc^ 

Decay of a single dressed state atom transforms an 
atomic symmetric Dicke state \N;Nq) into a combina¬ 
tion of the original state \N;Ne)^ a symmetric Dicke 
state \N;Ne — 1) with one fewer excitation, and N 
different other Dicke states {\N — l;Ne — l):-\l)i) in 
which the i-th atom is transferred to a non-qubit 
state (the qubit is lost), but which are still symmet¬ 
ric Dicke states for the remaining atoms. The result- 
ing state is y^\N; Ne)\0)p + ^/PdeNe\N; Ne - l)|l)p + 

V ^ \N - 1-, Ne - l)i\l)i\l)p where Pk = Xke~^’’ is 

the probability of losing/de-exciting {k = l/de) an atom 
over the cat creation time, with = 7(/c)y'^c (note that 
Ne ^ Y since the cat creation happens on the equator of 
Bloch sphere) and Pq = 1 — Pi — Pde- Here we focus on 
the regime where the probability of a single atom decay¬ 
ing is sufficiently small that the probability of two atoms 
decaying can be ignored. 

Tracing over the lost qubit and the photonic state one 

N 

obtains the density matrix pc = PopoP ^ X] P\P Pdepde^ 

i=l 

where po and pde are in the symmetric subspace with to¬ 
tal spin (J = y)^ while the p\ are in N different sym¬ 
metric subspaces with total spin (J = '^^^)- The po 
component corresponds to the ideal cat state. All the 
other components have very small fidelity with ideal cat 
states, primarily because the decay happens at a random 
point in time, which leads to dephasing. For example, 
de-excitation of an atom at {tde ^ [0,Tc]), leads to 

N 

IV-f (7e)) = 2-^/^ Y. VC{N,Ne) (15) 

Ne = l 

where F^^TVe-i) represents the dressed state energy of 
{Ne — 1) excited atoms, see Eq. Inserting the ex¬ 
pressions for EiVe and EjVe-i, one sees that de-excitation 
adds a linear term (iNeXotde) to the phase. This cre¬ 
ates a rotation around the 2 : axis on the Bloch sphere. 
The uncertainty in the time of decay tde therefore de¬ 
phases the cat state, resulting in the formation of a ring 
on the equator of the Bloch sphere, which has a small 
overlap with the ideal cat state. The fidelity of the re¬ 


sulting density matrix compared to an ideal cat state 
in the same subspace (which corresponds to the case 
where de-excitation happens at tde = 0) can be written 
as Pde = |(V’f®(7e)|'*/’c®(^rfe = 0))\'^dtde- When the 

size of the cat state is increased from = 10 to = 160, 
the fidelity of the generated cat in the de-excited sub¬ 
space is reduced from Fde = 0.2 to Pde = 0.045. The 
fidelity in each of the N subspaces where one atom was 
lost can be calculated in a similar way, yielding equiva¬ 
lent results. The total fidelity in the presence of Rydberg 
decoherence is then Fdc = Po P PlPl P PdePde ~ Tq. 

About 5% of Rydberg decoherence will transfer back 
to the excited state, which acts as dephasing (modeled 
by a Lindblad operator |e)(e|). The dephasing opera¬ 
tor commutes with the Hamiltonian for cat state cre¬ 
ation. Its effect can therefore be studied by having it 
act on the final cat state. For example, it can cause 
a sign flip of |e) for the first atom, resulting in a state 

I f^ \e)^g) ^jf^ \e)-^g) ^p_ 

plying the 7r/2 rotation results in a new energy cat 
-, which is clearly still a large superpo¬ 
sition in energy. So the effect of dephasing errors is rela¬ 
tively benign. Moreover, given the small relative rate of 
dephasing compared to loss and de-excitation, the prob¬ 
ability of having a sign flip over the cat creation time for 
the case with decoherence fidelity of Fde = 0.8 (consid¬ 
ered in Fig. will only be 1%. 

In conclusion, the fidelity of the cat state is, to a good 
approximation, equal to the probability of not losing or 
de-exciting any qubits over the cat creation time, Fde = 
Pq = g-0.95f 7erc_ 


V. ESTIMATE OF REALIZABLE CAT SIZE 

Taking into account the mentioned imperfections, 
Fig.i shows the achievable cat size as a function of the 
principal number n. Up to n ^ 80, the size increases 
with n. Higher n leads to a stronger interaction, hence 
allowing weaker dressing, and to smaller loss, favoring 
the creation of larger cats. However, for n ^ 80 the 
diminishing spacing between neighboring Rydberg levels 
(which scales like n~^) limits the detuning and hence the 
interaction strength, since xo = 2re^A and w has to be 
kept small, see Fig.[^ As a consequence, larger cat states 
cannot be achieved at higher principal numbers. 

Here we justify the behavior of Fig. in a more detailed 
scaling argument. For a constant fidelity the maximum 
achievable cat size N at each principal number n is lim¬ 
ited by Rydberg decay, Fde = where A = 0.95yTcXe. 
Let us analyze how A scales with N and n. The Rydberg 
decay rate scales as x\e) ^ w^n~^. In order to have a 
constant non-linearity fidelity of F^i =0.8, the dressing 
strength w has to scale like A-O-84^ ^he cat 

creation time Tc = ^ oc scales differently be¬ 

fore and after the transition point n ^ 80. Before the 
transition point the scaling of A can be obtained by not- 
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FIG. 4. (Color online) Maximum achievable cat size as a 
function of the principal number n of the Rydberg state. Ry¬ 
dberg state decay is adjusted to cause 20% infidelity. The 
interaction inhomogeneity is set to create less than 1% infi¬ 
delity, see Fig. 3, and the higher-order nonlinearities are set 
to create 10% (red circle), 20% (purple plus) and 30% (blue 
square) infidelity, see Fig. 2. The inset shows the required 
cat creation time as a function of n for the case where the 
higher-order nonlinearities cause 10% infidelity. 


ing that the trap size is a fraction of the blockade radius, 
A = oc where the exact value of the fraction 
coefficient is determined by Fjh^ see Fig. Therefore 
we conclude that A oc which states that before the 


transition point larger cat states are realizable by dress¬ 
ing to higher principal numbers, N (x for constant 
fidelity. However, after the transition point the small 
level spacing imposes a limit on the detuning, A cx n~^. 
Therefore after the transition point A cx which is 

independent of n. This prevents the realization of larger 
cat states at higher principal numbers. 

One sees that superposition states of over 100 atoms 
are achievable with good fidelity. In Fig. 4 the interaction 
inhomogeneity is tuned to create less than 1 % infidelity. 
Dressing to an S orbital is desired due to its isotropic 
interaction in the presence of trap fields. In Fig. after 
the transition point in n the detuning is chosen such that 
90% of the Rydberg component of the dressed state is 
bsns ^Si. Note that without a cryogenic environment the 
maximum achievable cat size in Fig. would be reduced 
from 165 to 120 atoms, see Sec. |VIIC[ 


VI. EXPERIMENTAL REALIZATION 


Experimental implementation of our scheme seems fea¬ 
sible. Rydberg excitations in Strontium have been real¬ 
ized over a wide range up to n = 500 [431446] . Rydberg 
dressing of two atoms has been used to create Bell-state 
entangled atoms m- Recently Rydberg dressing of up 
to 200 atoms in an optical-lattice has been reported [48] , 
where the collective interaction was probed using inter¬ 
ferometric techniques. Ref. [48| also identified a col¬ 
lective many-body decay process, which is however not a 


limiting factor for our scheme, as discussed in Sec. |VIID 

The Rydberg state bsns^Si is accessible from the 
bsbp^Po level with a 317nm laser field. The required Ry¬ 
dberg transition Rabi frequency fl^/27r (up to 15 MHz) 
can be obtained with a tunable single-frequency solid 
state laser of Ref. [49|. The relatively large detuning 
values (4MHz< A/27r < 340MHz in Fig. 4) make the 
interaction stable against Doppler shifts. 

Fulfilling the adiabaticity condition discussed in 
Sec. Inl is not difficult. In a highly adiabatic exam- 
ple, ^ = 0 . 01 , the dressing laser can be switched from 

zero to ^ = 15 MHz over 18 ns (for ^ = 270 MHz and 
165 atoms). For this example, 99.991% of the population 
returns to the ground state at the end of dressing, so adi¬ 
abaticity is almost perfect. This adiabatic switching time 
of 18 ns is many orders of magnitude shorter than the re¬ 
lated cat creation time of 1.4 ms. Adequate interaction 
timing precision is also required to align the created cat 
on the equator of Bloch sphere as explained in Sec. |H C 


For the 165-atom cat state mentioned above, a timing 
of order 5tc = ^ ~ 7.5 ns is re- 


precision 




quired for a phase precision of order 5(j) = = 7r/150. 

The Husimi Q function can be reconstructed based on 
tomography, i.e. counting atomic populations after ap¬ 
propriate rotations on the Bloch sphere. Modern fluores¬ 
cence methods can count atom numbers in the required 
range with single-atom accuracy [IHlIlo]. 


VII. OTHER SOURCES OF IMPERFECTION 


A. Effects of Atomic Motion in the Optical Lattice 


Laser manipulation of the atomic state leads to phases 
that depend on the atomic position. Atomic motion 
could therefore lead to decoherence. To suppress this 
effect, in the present proposal the atoms are confined to 
the ground states of the optical lattice traps. As a conse¬ 
quence, all position-dependent phase factors are constant 
over the course of the experiment and can be absorbed 
into the definition of the excited states. We now explain 
these points in more detail. Let us consider the atom, 
and let us assume that it is initially in the ground state 
(zero-phonon state) of its optical lattice site. We will de¬ 
note the corresponding state |^)j|0)j. Applying the part 
of the Hamiltonian that is due to the laser to this state 
gives We 

can rewrite the position operator Xj as the sum of the 
constant position of the site of the trap {x^j) plus 
a relative position operator = 5 (dt -b dj), where 

s = yj spread of the ground state wave 

function, ujtr is the trap frequency and (dj,dt) are the 
phononic annihilation-creation operators of the jth atom. 
In the Lamb-Dicke regime (77 = ^ ^ 1 ) one can expand 
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the exponential to get 

The phase factor is constant over the course of 

the experiment and can be absorbed into the definition 
of the atomic basis states by defining \e')j = 

The Hamiltonian describing the laser excitation can now 
be written in the new basis |^,0)j, |e',0)j, |e', l)j as: 

0 Vte \ ( |^,0)j \ 

0 0 |e',0), (17) 

0 ^tr j y 1^ 5 f)ji j 

Starting from the spin and motional ground state \g^ 0)^, 
the probability of populating the state corre¬ 

sponding to the creation of a phonon, will be negligible if 
CleV ^ With the parameters that we considered in 
our proposal ^ 1 kHz, r] = 0.1, ^ ^ 400 kHz) [31] 
the populatiou of |e',l)j will be eight orders of magui- 
tude smaller thau the populatiou iu the motioual grouud 
state. 


B. Effects of High Density 

The relatively small lattice spaciug of order 200um 
might raise coucerus about molecule formatiou aud level 
mixiug. At high atomic deusities there is auother po- 
teutial loss chauuel, Rydberg molecule formatiou m- 
Molecule formatiou ouly occurs wheu the attractive po- 
teutial due to Rydberg electrou-ueutral atom scatteriug 
moves the two biudiug atoms to a very small separa- 
tiou (of order 2um), where the biudiug euergy of the 
molecules cau iouize the Rydberg electrou aud form a 
Sr^ molecule [52]. Without the mass trausport, stepwise 
decay or iouizatiou of the Rydberg atom is ruled out by 
the quautizatiou of Rydberg state, as has beeu discussed 
aud experimeutally tested iu m, because eveu at high 
deusities the small molecular biudiug euergy of uearby 
atoms is orders of maguitude smaller thau the closest 
Rydberg levels for all the priucipal uumbers. The occur- 
reuce of iou pair formatiou is also highly uulikely iu this 
system [52]. We propose that coufiuiug the atoms by au 
optical lattice cau preveut the described mass trausport 
aud completely close the molecule formatiou loss chau¬ 
uel. High atomic deusity cau also lead to stroug level 
mixiug at short distauces [SSlISll- However, the experi- 
meut of Ref. m shows that the plateau-type iuteractiou 
cau persist iu the preseuce of stroug level mixiug because 
most molecular resouauces are ouly weakly coupled to 
the Rydberg excitatiou laser. 

C. Effects of Blackbody Radiation 

Blackbody radiation (BBR) could reduce the lifetime 
by transferring the Rydberg state population to neigh- 



FIG. 5. (Color online) Depopulation of Strontium Ryd¬ 
berg levels due to blackbody radiation (BBR) induced tran¬ 
sitions. a) BBR-induced transition rates (Einstein-B coef¬ 
ficients) from 5s80s^Si to the neighboring 5snp^P2 (dark 
blue), 5snp^Pi (light blue), 5snp^Po (blue) levels. The sum 
of these transition rates gives the total BBR-induced depop¬ 
ulation rate Tbbr- The inset is a 20 times enlarged view, 
b) Rydberg depopulation rates due to spontaneous decay (Tg 
shown in blue diamond) and BBR-induced transitions (E bbr) 
at environment temperatures of 300K (red circle), 95K [56] 
(purple circle), and 3K [53 (green circle) as a function of the 
principal number. The use of a cryogenic environment signif¬ 
icantly suppresses the unwanted effects of BBR. 


boring Rydberg levels (with both higher and lower prin¬ 
cipal numbers n) as illustrated in Fig. The BBR- 
induced transition probability is given by the Einstein B- 
coefBcient Tbbr = nP — [SlHin], where 

/ / e -1 

T is the environment temperature, is the Boltzmann 
constant and both ujif and Aif are defined in Sec. |IV Aj 

At the environment temperatures of 300K, 95K |56j 
and 3K [57], including the BBR-induced transitions in¬ 
creases the total decoherence rate Fg by 120%, 40% and 
1% (see Fig. If) for n ~ 80, which results in maximum 
achievable cat sizes of 120, 150 and 165 atoms respec¬ 
tively (considering F^i = 0.7, = 0.8). Note that 

cryogenic environments with 95K and IK were used in 
a Strontium lattice clock experiment [56] and in a cavity 
QED experiment with Rydberg atoms m respectively. 

BBR could also disturb the Ramsey-type interferom¬ 
etry used for detecting energy decoherence by produc¬ 
ing an AC stark shift; this effect is quantified in section 
|vmi Furthermore, BBR-induced decoherence could be 
inhomogeneous due to temperature inhomogeneities in 
the environment. This would introduce unwanted cou¬ 
pling to non-symmetric Dicke states in the cat creation 
process. The use of a cryogenic environment significantly 
suppresses these effects as well. 


D. Effects of Collective Many-body Decoherence 

BBR-induced transitions to neighboring Rydberg lev¬ 
els (see Fig. |5^ can also lead to collective many-body 
decoherence |48l [58] . The interaction between the tar- 



























get nS Rydberg level and some of the populated neigh¬ 
boring n'P levels is of a strong long-range dipole-dipole 
type due to the formation of Forster resonances. This 
strong interaction causes an anomalous broadening [58] . 
The mentioned decoherence process only starts after the 
first BBR-induced transition occurs. However, the weak 
dressing strength and small ensemble size {N < 200) in 
our scheme make the probability of populating the tar¬ 
get Rydberg state and consequently neighboring Ryd¬ 
berg levels very small. For example at the environment 
temperatures of 300K, 95K and 3K and for dressing to 
n ~ 80, the probabilities of not populating the strongly 
interacting neighboring Rydberg levels over the cat cre¬ 
ation time for cat sizes of 120, 150 and 165 atoms re¬ 
spectively are Pbbr{^) = exp{—Y'i^‘^^BBR^c) =98.63%, 
99.26% and 99.96% respectively. It has been observed 
in the realization of many particle Rydberg dressing [48] 
that when the transition probability is low enough (of 
the order of Pbbr{0) ^ 82%, as can be calculated from 
the information provided in Ref. [48]) the many-body 
decoherence effects are negligible and decoherence rate 
is dominated by the Rydberg depopulation rate (see 
Sec.|iV). 


VIII. TESTING ENERGY DECOHERENCE 

In the context of modifications of quantum physics, 
decoherence in the energy basis is quite a natural pos¬ 
sibility to consider IMig. It is usually introduced as 
an additional term in the time evolution for the den¬ 
sity matrix that is quadratic in the Hamiltonian, ^ = 
^ [H, p\ — [i^, [i^, p]], which leads to a decay of the off- 

diagonal terms of the density matrix in the energy basis 
according to pnm{t) = [18], where 

'jE = Here uJnm IS related to the energy differ¬ 

ence of the two componants and a can be interpreted as 
a timescale on which time is effectively discretized, e.g. 
related to quantum gravity effects. It is of interest to es¬ 
tablish experimental bounds on the size of a, which could 
in principle be as small as the Planck time (10“^^ s). 

The corresponding decoherence rate for the energy cat 
in this proposal would he je = <t( ^^^ )^, where AE 
is the energy difference between the ground and excited 
state of each qubit, and N is the cat size. To detect the 
energy decoherence one prepares the energy cat state, 
followed by a waiting period. To observe the decoherence 
effect, one detects the Ramsey fringes for the revival. 
The visibility of the Ramsey interference is also sensitive 
to other decoherence sources, where in the absence of 
dressing laser the dominant ones are the trap loss rate F, 
which reduces the visibility by a factor exp(—VFt), and 
phase diffusion that is explained below. 

The large energy difference of the cat state increases 
the sensitivity of the Ramsey interferometry that we are 
using for the detection of energy decoherence. There¬ 
fore, it is important to consider the effect of fluctuations 
in the detuning between the laser and the atomic tran¬ 


sition. Let us first note that the cat state is more sensi¬ 
tive to multi-particle (correlated) than to single-particle 
(uncorrelated) noise, which results in a phase diffusion 
affecting the visibility of Ramsey fringes by and 

respectively [59]. Comparing the two cases, 
correlated fluctuations should be ^/N times more sta¬ 
bilized than uncorrelated fluctuations. The most impor¬ 
tant source of noise in our system is the fluctuation of 
the laser frequency. A probe laser linewidth as narrow 
as 26 mHz m has been achieved in optical atomic clock 
experiments, and there are proposals for much smaller 
linewidths HU Eg with recent experimental progress [63] , 
justifying our example of a lOmHz linewidth, see below. 


Other sources of multi-particle and single-particle noise 
have been well studied in the context of Strontium atomic 
clocks [64| [65] and are comparatively negligible. Here 
we address a few of them in our scheme. One of the 
noise sources is the trap fields intensity fluctuation; how¬ 
ever, using the magic wavelength makes the atomic tran¬ 
sition frequency independent of the trap laser intensity. 
Considering the variation of the Stark shifts due to the 
trap laser as a function of frequency at the magic wave¬ 
length [66], the relative scalar light shifts could be kept 
within O.lmHz uncertainty by applying a trap laser with 
a IMHz linewidth. In addition to the scalar light shift, 
the inhomogeneous polarization of trap fields in 3D op¬ 
tical lattices can result in an inhomogeneous tensor light 
shift [67]; however, the use of the bosonic isotope ^^Sr 
with zero magnetic moment cancels the tensor light shift 
[68] in our scheme. Environmental temperature fluctua¬ 
tions (ST) also lead to atomic frequency fluctuations that 
are proportional to T^5T due to the BBR-induced light 
shift [64]. This is another reason why a cryogenic en¬ 
vironment is advantageous. For example controlling the 
environment temperature of 95K [56] to within a range 
of ST = IK keeps the BBR-induced noise shift below I 
mHz. 


A conservative estimate of the experimentally measur¬ 
able energy decoherence rate can be obtained by consid¬ 
ering the case where the energy decoherence dominates 
all other decoherence sources during the waiting period. 
Increasing the cat size N is helpful because it allows one 
to enhance the relative size of the energy decoherence 
contribution. For example, choosing t oc keeps the 
loss and phase diffusion contributions fixed, while the 
energy decoherence still increases proportionally to N. 
Using a cat state with N = 165 atoms (see Fig. 4), 
which corresponds to NAE = 300 eV, assuming a laser 
linewidth of 10 mHz (see above), and considering a trap 
loss rate of F = 10 mHz m, the minimum detectable dis¬ 
cretization time scale a is of order 10“^^ s. This would 
improve the measurement precision by 4 and II orders 
of magnitude compared to what is possible based on Ref. 
[26] and Ref. [25] respectively. 
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Appendix A: Effects of Interaction Inhomogeneity 


Here we explain the steps in calculating the effects 
of inhomogeneous interaction on the cat st ate’s f idelity 
Fih = |(^ = I 77 = l)p (see Sec. IIIB). Tay¬ 
lor expanding and considering the expectation 

values (7? = = 1) = 1/2 and (7? = = 


1 ) = \ one obtains an estimate for the fidelity. 

The first order of the expansion is zero because we de¬ 
fined €ij as fluctuations around a mean value. The sec¬ 


ond order can be calculated using {r] = l\V^fj\r] = 1 ) = 
I Cijim^ijeim, where Cijim = 1/16 if all the in- 

i^j l^m 


dices are unequal, Cijim = 1/8 if there is a pair of equal 
indices, and Cijim = 1/4 when there are two pairs of 
equal indices. The convergence of the expansion for the 
fidelity can be tested numerically. In Fig. 3b of the paper 
the ratio of the third order to the second order of the ex- 
pansionforFm is = 10 “®,5xl0-^SxlQ-^,SxlQ-^ 

for ^ = 0.1,0.2,0.3,0.4 respectively, suggesting good 
convergence in this regime. 


Appendix B: Dipole Matrix Elements 

In Strontium one needs to consider both valence elec¬ 
trons (|i) = \niin 2 iliil 2 iLiSiJiMi)) in the calculation of 


the dipole matrix elements m 


|(i|r|/)p = max 

j/ 1 

Li S Lf 


{hi, hf) {“^Lf + 1)(2 J/ + 1 )( 2 L, + 1 ) 

{ hi hf } \^^^Mr\n2fhf)^,^^ 


where L and S are the total orbital angular momentum 
and spin, I and s refer to individual electrons, and J and 
M refer to total angular momentum. The active electron 
in the transition is labeled by 2 , and (n2i/2ikl^2/^2/) is 
the radial dipole matrix element between initial and final 
state, and the curly bracket is a Wigner- 6 j symbol. 
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